clear;clc
% x0=pi/4;x1=pi/2;
% k=0;counter=0;
% while k==0
%     counter=counter+1;
%     x=(x0+x1)/2;
%     y1=4*cos(x);
%     y2=exp(1)^x;
%     statu=y1-y2;
%     if statu>0
%         x0=x;
%     else
%         x1=x;
%     end
%     if abs(statu)<0.001
%         k=1;
%     end
% end

% f=@(x)exp(1)^x-4*cos(x);
% ff=@(x)exp(1)^x+4*sin(x);
% x=pi/4;
% k=0;
% counter=0;
% while k==0
%     counter=counter+1;
%     x0=x-f(x)/ff(x);
%     if abs(f(x0))<0.001
%         k=1;
%     end
%     x=x0;
% end

f=@(x)exp(1)^x-4*cos(x);
x0=pi/4;x1=pi/2;
k=0;
counter=0;
while k==0
    counter=counter+1;
    x2=x1-f(x1)*(x1-x0)/(f(x1)-f(x0))
    if abs(f(x2))<0.001
        k=1;
    end
    x1=x2;
end

% f=@(x)exp(1)^x-4*cos(x);
% options=optimset('TolX',0.001);
% fzero(f,[pi/4,pi/2],options)



